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The inclusion of charging and spin-exchange interactions within the Universal Hamiltonian description of 
quantum dots is challenging as it leads to a non-Abelian action. Here we present an exact analytical solution 
of the probem, in particular, in the vicinity of the Stoner instabilty point. We calculate several observables, 
including the tunneling density of states (TDOS) and the spin susceptibility. Near the instability point the 
TDOS exhibits a non-monotonous behavior as function of the tunneling energy, even at temperatures higher 
than the exchange energy. Our approach is generalizable to a broad set of observables, including the a.c. 
susceptibility and the absorption spectrum for anisotropic spin interaction. Our results could be tested in 
nearly ferromagnetic materials. 
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The physics of quantum dots (QDs) is a focal point 
of research in nanoelectronics. The introduction of the 
"Universal Hamiltonian" [Tl[2] made it possible to sim- 
pUfy in a controlled way the intricate electron-electron 
interactions within a QD. This provided one with a 
convenient framework to calculate physical observables. 
Within this scheme interactions are represented as the 
sum of three spatially independent terms: charging, 
spin-exchange, and Cooper channel. Notably, even the 
inclusion of the first two terms turned out to be non- 
trivial: the resulting action is non-Abelian [Sill]. At- 
tempts to account for those interactions in transport 
involved a rate equation analysis [SJIB] and a perturba- 
tion expansion [4] . Alhassid and Rupp [5] have analyzed 
some aspects of the problem (see below) exactly. It is 
known that in the presence of significant spin-exchange 
interaction such systems can become Stoner unstable. 
More precisely, one distinguishes 3 regimes of behavior 
as function of increasing the strength of the exchange in- 
teraction: paramagnetic (no zero field magnetization), 
mesoscopic Stoner regime (finite magnetization whose 
value increases stepwise with the exchange) and ther- 
modynamic ferromagnetic phase (magnetization is pro- 
portional to the volume) [2]. Both the mesoscopic and 
thermodynamic phases manifest (Stoner) instabilities 
towards ferromagnetic ordering. The presence of en- 
hanced quantum and statistical fluctuations underlying 
such instabilities calls for a full-fledged quantum me- 
chanical treatment of the problem. 



Here we present an exact analytic algorithm to tackle 
this challenging problem. We employ our approach to 
a few physical variables within the mesoscopic Stoner 
regime, but it can be used to tackle the broad range of 
problems involving spin and charge on a QD, and be 
extended to the thermodynamic ferromagnetic regime 
too. As examples we calculate the following quanti- 
ties: the partition function, the magnetic susceptibility, 
the distribution function of the total spin, the tunnel- 
ing density of states (TDOS), and the sequential tun- 
neling conductance. Our approach allows us to obtain 
analytic results as one approaches the Stoner instabil- 
ity. Below we list possible applications of our method 
to other physical observables and extensions beyond the 
Universal Hamiltonian. The physics discussed here can 
be best tested in quantum dots with materials which 
are close to the thermodynamic Stoner instability, e.g., 
Co impurities in Pd or Pt host, Fe dissolved in various 
transition metal alloys, Ni impurities in Pd host, and Co 
in Fe grains, as well as new nearly ferromagnetic rare 
earth materials [71[H]- 

The main reason why, in this context of a QD, the 
treatment of the exchange term is non-trivial, is the 
non-Abelian nature of the action. One needs to tackle 
time ordered integrals of the form 



Texp 



(1) 
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Here 6p is a dynamical, quantum held operating on the 
spin s-y (whose x component is proportional to the Pauli 
matrix ax etc.); p and 7 are indices to be elaborated 
below; T is a time ordering operation. Wei and Nor- 
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man [5] , addressing the problem of a quantum spin sub- 
ject to a prescribed classical time-dependent magnetic 
field, have elegantly shown that by preforming a non- 
linear transformation from 8p,6y, 6^ to a set of other 
variables (cf. Eq. ((T5|) ). Eq. ([1]) can be written as a 
product of 3 Abelian terms (cf. Eq. ([T6))). Even so, that 
problem could not be solved. The problem of a quan- 
tum field appears to be even more intricate. To solve 
it we employ here a generalized Wei-Norman-Kolokolov 
(WNK) method [TO]. 

We consider a quantum dot of linear size L in the 
so-called metallic regime, whose dimensionlcss conduc- 
tance gxh ~ ETh/S ^ 1. Here Etyi is the Thoulcss 
energy and S is the (spinless) mean single particle level 
spacing. We account for the following terms of the Uni- 
versal Hamiltonian 



H = Ho + Hc 



(2) 



Here, denotes the spin (a) degenerate single particle 
levels. The charging interaction He ~ Ec{n — Nq)'^ ac- 
counts for the Coulomb blockade, with n = = 



y 



being the particle number operator; A^o 



represents the positive background charge. The term 
Hs =■ —JS^ represents spin interactions within the dot 



{S = Y^a^a = ^Y.a^i,cr'^'^'^'"-a,T'), with the Compo- 
nents of cr comprising of the Pauli matrices. 

The imaginary time action for this system reads: 



dr. 



Here fi is the chemical potential, /3 = 1/T, T the tem- 
perature, and we have introduced the Grassmann vari- 
ables *a = {i'a^.i^aiV =■ (V'aT'^a;) ^ rCprCSCnt 

electrons on the dot. 

Employing a Hubbard-Stratonovich transformation 
leads to a bosonized form 



5r — Co + M + if/) + 



cr $ 



$2 02 
+ — iNQ(h 



where and $ arc scalar and vector bosonic fields re- 
spectively. The SU{2) non- Abelian character of the ac- 
tion poses a serious difficulty. It prevents one from per- 
forming a gauge transformation which works effi- 
ciently in the Abelian U (1) (charging only) case pTHTO] . 
Employing the Wei-Norman-Kolokolov trick we are able 
to overcome this difficulty. 



Results. — Below we present our main results. The 
TDOS is given by the following exact expression 



1 + 

Z 



■ E 



^-l3E^{n-Noy+l3tm+l3Jm{m+l) 



E - + M - £^c(2n - 2Ao + 1) - J{m + 

a 



X < 2m 



+(2m + l) Z„,Z„^( 



(4) 



Here n^^{n:^) represents the number of spin- up (spin- 
down) electrons, the total number of electrons n = 
+ ri_i_, m = (n^^ — n^jl. Note that for m ^ (m < 0) 
the total spin S = m (S = ~m — 1) respectively. The 
normalization factor 



Z =^(2m- 



^)Zn^Zn^e 



-l3[Ea{n-No) -tJ.n-Jm{7n+l}] 



(5) 

coincides with the grand canonical partition function 
for the Hamiltonian ^ [5]. The quantity Z^ = 



/o 



2ii' d0 ^-iON 



2-K 



(l -I- e*^ '^'^t) is the canonical partition 



function of N noninteracting spinless electrons, and 



^w(ea) = /o 



— f2T dS „~idN 



n 



1- 



) dete 



rmmes 



the canonical partition function of a system of A'^ non- 
interacting spinless electrons under the constraint that 
level a is not occupied. 

Eqs. (jH) and ([S]) allow us to study a host of physical 
observables for a given spectrum of single-particle levels 
{eq}. At low temperatures, T < these observables 
are sensitive to details of the spectrum; their statistical 
averages would depend on the symmetry group of the 
spectral distribution [14j . 

We now discuss a few quantities of interest. The 
static spin susceptibility can be computed as x = 
(l/3)(91nZ/(9J. At high temperatures, 5 < T < 
^/ ln( Ji,/T), Ji, — J5/{5 — J), the average static spin 
susceptibility is given by 



1 1 



1 



X = 77 



25- J 12T{5- Jf 



This expression, underlining the divergence at the 
Stoner instability point, differs from that found by Kur- 
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land et al. [5] and by Schechter ^ ^\ Near the 
Stoner instability, S — J S, it is the first (second) 
term of Eq. (O that dominates when T J^, {T <^ J^,) . 
For T ^ Ji, the susceptibility behaves like a param- 
agnetic Fermi liquid (with an upward renormalized g- 
factor). As the system is driven towards the Stoner 
instability limit one crosses over to the low temperature 
regime, T <C J*, and a non-Fermi liquid (Curie) behav- 
ior, sets in, x ~ {S'^)/T, where the average spin scales as 
a/ (5^) ^ J*/<5. Note that the latter approximates the 
discontinuous growth of the ground state spin of a spe- 
cific single electron spectrum (e.g. uniformly spaced), 
when J/ 6 is increased in the mesoscopic Stoner regime 
towards 1. No dynamical spin response x(w 7^ 0) exists 
unless the dot is connected to reservoirs, or anisotropic 
spin interaction is considered. 

The average moments of the total spin can be 
found from the partition function Z as ([5'^]'^) — 
T'^Z^^d'^Z/dJ^. It can be characterized by the dis- 
tribution function of S^, Vs^{x) which can be found 
from Eq. ([5]). Near the Stoner instability 6 ~ J 6, 
and for the same range of temperatures as in Eq. ([5]), 
the distribution becomes 



v(e)/V(, 



Vs^ix) = 2J^e-'3■^'/4sinh(/3,5^/J) e 



(8) 



The broad asymmetric non-Gaussian nature of the dis- 
tribution becomes manifest in the high temperature 
limit, and is not due to statistical fluctuations of the 
single particle levels but rather due to the effect of the 
exchange interaction. 



^^The average static spin susceptibility has been calculated in 
Ref. [2] near the Stoner instability, 5 — J <g 5. In our notations, 
the result of Ref. [2] at T > J* becomes (sec Eqs.{4.8), (4.13b), 
(4.15) of Ref. [2]) 



CO 



J 



1 ■ 



- Cl—= +C2 — 
y/T T 



where numerical coefficients cq = 1/3, c\ = y/iv/A, and C2 ~ 0.238 
for unitary ensemble and cq = 1/3, ci = \/27r/4, and C2 ~ 0.227 
for orthogonal ensemble. The result of Ref. [2] contradicts our re- 
sult ||6jl in which cq = 1/2, ci = and C2 = 1/6 are independent 
of the ensemble statistics of the single-particle levels. According 
to Ref. [2], at T = 0, (see Eq.(4.19) of Ref. [2]) x [5/(<5 - J)]^. 
As one can see from Eq. l(6]l, our result for T <^ J* smoothly 
interpolates into the result of Ref. [2] for T = 0. 

'''Our result for x implies that the magnetic field tends to 
zero first (before, e.g., temperature). The result found by 
Schechter 1151 is valid in the limit of vanishing temperature but 
finite magnetic field (provided an additional coarse graining is 
performed). Generalization of Eq. ^ to finite magnetic field 
resembles the result of Schechter at magnetic fields larger than 
temperature | 14| . 
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Fig. [U TDOS in the Coulomb valley. The solid 
(dashed) line corresponds to J/(5 = 0.92, (5/T = 0.35, 
and jJjT = 3.95 (J/(5 = 0.92, 5/r = 0.95, and 
J*/r = 10.70). The inset depicts the nonmonotonic 
behavior. 

We next consider the average TDOS at (5 < T. The 
most interesting regime seems to be that of interme- 
diate temperatures, T <C J*. Under the assumption 
/J, ^ Tin J,t/T, Eq. ^ can be simplified, leading to 



E(£) 



n,(7— it 



2j/l J. '^-^^ 




,-f)E^{n-Nor 



(9) 



Here Q'^ = Edji -No + cr/2), is the averaged TDOS 
for noninteracting electrons, and 



1 / 7ra;\ _v(^_if+v^cos'^. 



•^(2;, y)= 2 '^s" v*^°^ "2" ' ^ 



(10) 



1 - $ 



,^(x-\x\) 



X XI (-l)"e-^l^l'"+^™("+i)6'(|x| - 2m - 1). 

6{x) is the Heaviside step function (0(0) = 0), and the 
error function $(z) = (2/y^) /J^ exp(— t^)(it. As x is 
varied for a fixed y, J-{x, y) exhibits damped oscillations 
with a period 4 (equivalent to an energy scale 4J^). In 
the limit y ^ 1 considered here, these oscillations are 
strongly suppressed, and only the first maximum re- 
mains visible. It leads to the appearance of a maximum 
in the TDOS as illustrated in Figs. [1] and H The scal- 
ing of these oscillations with ^ {S"^) ^ J*/i5 indicates 
that they are due to precession of the spin of the in- 
jected electron about the effective magnetic moment in 
the dot. This additional structure in the TDOS reflects 
enhanced electron correlations due to the exchange in- 
teraction. At higher temperatures, T 3> J*, there is no 
interesting signature of spin-exchange on the TDOS. 
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Fig. [21 TDOS at the Coulomb peak. The parameters 
are the same as in Fig[T] The insets depict the non- 
monotonic behavior. 

One can compute the sequential con- 
ductance through the QD employing G = 
Go / de{—d f f{£) / de){V{e) / vq) , where Go is the 
conductance of the non-interacting QD. The maximal 
value of G will be enhanced by a factor 1 + J/2J* 
due to the exchange term. Much more interestingly, 
the non-linear conductance at the Coulomb peak will 
exhibit non-monotonic behavior, similar to Fig. [2] |14j . 

Derivation. - Below we describe the main steps of 
the derivation. Further details will be given in |14j . 

The TDOS, v{e) = -(l/7r)ImE„,, G«,(^), is de- 
termined via the imaginary part of the retarded Green's 
function, G^,{t,t') = -^^^(^ - i')({ao,a(i), at _,(i')}>ff 
of the Hamiltonian ([2|). The imaginary time 
Green function is given by Gaa{Ti,T2) = 

- (Tr -00 ,a (ti ) ,a (t2 ) ) Stot • 

The exact one-particle Green function for the Hamil- 
tonian (12) can be written as 



Gacr{Tl,T2} = J d(j)o 



z 



-D(ri2,0o)Sao-(ri2,0o), 



-ttT 



(11) 

where T12 — ti — r2, 00 is the static component 
of 0, the grand canonical partition function Z = 
JI'^t'^'^o D{0,(j)o)Z{(l)o), and the so-called Coulomb- 
boson propagator reads pUfT^ 



,i4>o{l3k+T)-fiE,Xk-N„f-2E^{k^No)T 



The one-particle Green function Qaa{'Ti,T2,4>o) ap- 
pearing in Eq. ([TT|) is defined as ^a(j(Ti, T2, (/)o) = 
— (^Ti^a.a{'Ti)i}}a.a{'T2))s- Avcragc is taken with respect 
to the action S ~ Jo J2a^ o'dr'^ a — H ■ Here 
H ~ Hq + Hs with Hq in which is replaced by 
Eq = Eq — A* + *0o- Remarkably, the charge and spin 



degrees of freedom are almost disentangled in the action 
S. The latter involves only the spin-interaction part of 
the Hamiltonian ([2|) . Traces of the charging- interaction 
are encoded in the variable (f>o, leading to a small imag- 
inary shift of the chemical potential. Subsequently, the 
one-particle Green function can be written as 

Gaain > T2) = -Z^^]Cacr{-iTl2, -«Ti2 + 

lC^a{t+,t-) = Tre-**+«al^,e'*-«a„,, (13) 

and Z{4io) = Trexp(— /SH). In order to evaluate 
the trace we perform Hubbard-Stratonovich transfor- 
mations of the terms e^'**^^ in the evolution operators 
and obtain 



/c,,(i+,t_) = H fv[ep] 



dt'el 



(14) 



X Tr 



Here is defined in Eq. 



(H]). We have defined 
the bosonic fields Op, p — ±. In order to employ the 
WNK trick we use a Hamiltonian evolution of our op- 
erators rather than a path untegral representation of 
Q. Note that while H is time independent, the factors 
A^'^ involve time ordering (T). This is due to the non- 
commutativity of the spin-operators Sj. 

In order to overcome the intricacy of time-ordering 
we use the following transformation of variables I16| in 
the functional integral in Eq. [TU] . 



Pp 2K^Kp 



-ipkP + ppKl - {nP^Kp P, (15) 



which recasts the time-ordered exponent as a product 
of simple Abelian ones: 



(16) 



X exp 



dt'KpP{t')e-'PJ'o '^^Pp(^Ut' 



Here we employ the initial condition k^(0) =0 0, and 
~ ± is^. We stress that Eqs. (fT^ and are 
valid for a general spin operator. In order to preserve the 
number of field variables (three) we impose the follow- 
ing constraints on the otherwise arbitrary new complex 
variables: pp = —p*p and k+ = ('^p )*• The quantity 
)(^aa{t+,t-) can be then evaluated as 



n 



V[pp,KP]e' 



xeT.o --p-)C(t+,t_)[]6^(t+,L), (17) 
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with Caa and B-y given in terms of single-particle traces: 



r — tr 



B-y = tr 



A^-\t^) . (18) 



The expression for Z can be obtained from Eq. (|17|) by 
the substitution of Ba for Caa- We can now evaluate 
the single- particle traces in Bj and Caa- The fields 
appear in Bj. It turns out that the integration over 
first, and then pp, can be performed exactly, yielding 



IC. 



e t 



■at 



dh 



sinh(/3/i) n J] [ 



1 + e 



X 2J e*^°*= + ^-e J:^^? {2ph + isJt^s)- (19) 

Next, we perform the integration over h in Eq. (|19p. 
substitute it into Eq. ((T3)) and calculate the exchange- 
only Green function, Qaa- Then, integrating over (pQ in 
Eq. (fTT|) we obtain the full Green's function Gaa- Em- 
ploying the general expression [18j 



- — cosh — 



at 



it- 



, (20) 



we, finally, find the TDOS (U). In a similar way we 
obtain the partition function Z ([5|). 

Within WNK method one may still have some free- 
dom in selecting regularization of the functional inte- 
grals. It is thus useful to check the validity of our re- 
sults against some benchmarks. Our non-trivial checks 
are: i) Eq.(5) for Z agrees with the exact derivation 
in Ref. [S]. ii) The TDOS g]) satisfies the sum rule: 
/ dev{e)fF{e) = Td\nZ/dfj. HI], iii) For J = our 
results for the TDOS coincide with those of Ref. (T5] . 
iv) Our results for Z and z/(e) agree with a direct cal- 
culation for single and double level QDs. 

In summary, we have addressed here the interplay of 
charging and spin-exchange interactions of electrons in 
a metallic quantum dot. Even within the simple Univer- 
sal Hamiltonian framework, this problem leads to a non- 
Abelian action, and necessarily requires the evaluation 
of non-trivial time-ordered integrals. Our method is ap- 
plicable to the vicinity of the Stoner instability (well in- 
side the mesoscopic Stoner unstable regime), and could 
be extended to the ferromagnetic regime. Other ex- 
tensions include the study of anisotropic spin-exchange 
(where the non- vanishing a.c. susceptibility, absorption 
and TDOS are of particular interest), cotunneling con- 
ductance, and an explicit inclusion of the leads. 



As a demonstration of the usefulness of our exact 
solution we have calculated several quantities: the par- 
tition function, the magnetic susceptibility, the distri- 
bution function of the spin, the TDOS, and the linear 
and non-linear conductance at the Coulomb peak. Some 
of these quantities are amenable to experimental tests. 
Examples: the broad distribution of the spin would im- 
ply significant sample-to-sample fluctuations of the mea- 
sured susceptibility; the latter can be used to determine 
the distance (1 — J/ 5) from the Stoner instability; the 
relative magnitude of the predicted non-monotonicities 
in the TDOS and the conductance may exceed 5 — 10% 
in materials close to the Stoner instability such as Pd 
{J/S = 0.83) or YFe2Zn2o {J/S = 0.94) 0. 

Previously, Alhassid et al. have calculated exactly 
the partition function, matrix elements of al^^^aaa [5], 
and many-body eigenstates which are also eigenstates 
of the total spin operator [T^ . That approach could be 
employed for the calculation of other observables. Our 
independent approach is more manageable for the cal- 
culation of higher correlators, the inclusion of exchange 
anisotropy, as well as to further generalizations, as in- 
dicated above. 
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